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Abstract 

The scaling of the time delay near a "bottleneck" of a generic saddle-node bifurcation is well- 
known to be given by an inverse square-root law. We extend the analysis to several non-generic 
cases for smooth vector fields. We proceed to investigate C° vector fields and demonstrate that 
parameter nonlinearities and phase-space nonlinearities have distinct effects. Our main result is 
a new phenomenon, which we call "quantitative bifurcation" in two-parameter families having a 
saddle-node bifurcation upon changing the first parameter. We find distinct scalings for different 
values of the second parameter ranging from power laws with exponents in (0, f ) to scalings given 
by O(l). We illustrate the quantitative bifurcation by a an overdamped pendulum with varying 
length. 
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Saddle-node bifurcations have been extensively studied in dynamical systems. The normal form in 
the context of ordinary differential equations is 

x = r + x 2 for x £ K and rel (1) 

where r is the bifurcation parameter. Solving for the fixed points we set r + x 2 = 0, which has two 

solutions x± — ±i/^r giving no fixed points for r > and a non- hyperbolic fixed point x± = at 

r = 0. For r < we obtain an attracting fixed point x_ and a repelling fixed point at 2;+, hence a 

saddle-node bifurcation occurs at r = 0. If we set x(0) = xq we can solve ([T]) by separation of variables 

and obtain 

r x(t) ^ rt 

xds = dt (2) 

r + s z Jr. 



Using trigonometric substitution we get 
00 . 
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tan- 1 P^ - tan- 1 ^ = VF (3) 



Therefore the time a trajectory spends in the interval [—1,1] is 

= t (4) 



2tan- 1 (^ 



Taking r — > + gives 

t~4 ( 5 ) 

In the literature this scaling law is referred to as intermittency, "bottleneck" or saddle-node ghost. 
One of the earliest references is [7j where saddle-node bifurcations of maps are investigated. Textbook 



references include [3] and [5]. Applications of the scaling law to physical systems can, for example, be 
found in [1], [3] and [8]. We remark that the square-root scaling is found for any generic saddle-node 
bifurcation in sufficiently smooth vector fields independent of the dimension of the phase space. To 
justify this statement, recall the following theorem (see [3]): 

Theorem 1.1. Consider y — g(/i, y) with /ieR and y £ W l . Assume that for /i — fiQ there exists an 
equilibirum yo such that: 

1. D y g(fio,yo) has a simple eigenvalue with right eigenvector v and left eigenvector w. 

2. D y g(fiQ,yo) has k eigenvalues of negative real part and n — k — 1 eigenvalues with positive real 
part. 

3. w((dg/dfj,)(y ,Ho))^0 

4- w{D 2 y g{^,ya){v,v)) ^0 

then there is a smooth curve of equilibira in R x W. n passing through (/xq, J/o) tangent to the hyperplane 
{/xo} x W 1 with no equilibrium on one side of the hyperplane for each value \x and 2 equilibria on the 
other side of the hyperplane for each \x. The two equilibria are hyperbolic and have stable manifolds 
of dimensions k and k + 1 respectively. 

Furthermore the set of equations satisfying Theorem 11.11 is known to be open and dense in C°° 
one-parameter families with an equilibrium having a simple zero eigenvalue. 

If we assume in addition that D y g(po, yo) has no eigenvalues with positive real part we obtain that 
the dynamics near (/zo,yo) has n — 1 exponentially attracting directions outside of a 1-dimensional 
center manifold M . On M the dynamics is given up to a normal form transformation by x = /i + x 2 
for x G M (see e.g. [5]). 

Hence we can restrict to the equation 

x = f(r, x) for x € R and r € 1R 

We shall from now on consider one-dimensional flows only. Necessary conditions for a saddle-node at 
(r, x) = (0, 0) are /(0, 0) = = f x (0, 0). The generic sufficient conditions are given by: 

fxx (0,0) ^0 (6) 

fr (0,0) ^0 (7) 

We can easily see, using the implicit function theorem for f x at (0, 0), that f{x, 0) is locally quadratic. 
Furthermore / r (0, r) is locally linear, therefore the square-root scaling law is not affected by a change 
of coordinates transforming the equation x — f{r,x) into the normal form (TTJ). Hence all smooth 
vector fields with generic saddle-node bifurcations scale in the same way with respect to a bottleneck 
occuring for r — ► + . 

In the next section we briefly address the question, what happens to the scaling law when the 
saddle-node is non-generic. In this case either ((6]) or ([7|) are violated and we refer to such a situation 
as a degenerate saddle-node. 

2 Degenerate Saddle-Nodes 

We treat each of the cases for degenerate saddle-nodes in turn. 



2.1 Case 1 - / r (0, 0) = 

The linear crossing condition for r at r — fails and we consider the equation 

x = R(r) + x 2 for x £ R, r £ R and R £ C°°(R) 
without restrictions on the first derivative of R. In particular we can show: 

Proposition 2.1. Let a : (0, oo) — ► K + with a(r) = 1/r for k £ {2,3,...}. Then there exists a 
function R(r), R £ C°°(R), such that the equation x — R(r) + x 2 has a saddle-node bifurcation at 
r = 0. The scaling law for this saddle-node is given by t ~ a(r) as r — ► + . 

Proof. Define R(r) = l/a(r) 2 where we use the even extension for l/a(r) 2 to define R(r) on the entire 
real line so that R(r) = r 2k . Note that R(0) = and R £ C°°(R). Then fix any r £ R and set 
R(r ) — k > 0. Using the same separation of variables argument as in ©-flS) we consider 

x = k + x 2 

and obtain the time a trajectory spends in [— 1, 1] as 

2 tan" 1 (—') 
t = VV ^ = 2a(r ) tan" 1 (a(r Q )) 

so that since ro was arbitrary we get that t ~ a(r) = r fc as r — > + . D 

Remark: Note that we have included topolocially degenerate cases such as x = r 2 + x 2 in Propo- 
sition [2T] Disregarding these cases excludes all scaling laws with even powers k = 2to. 



2.2 Case 2- f xx (0,0) 
Consider the equation: 



F(x) for x £ R and r £ 



We assume for now that F £ C°°(R). Notice that we still require that F(0) = = F x (0) as necessary 
conditions for a saddle-node at (r, a;) = (0, 0). Also we assume without loss of generality that F(x) > 
for i^0. Considering the Taylor expansion of F at we get 

F(x) = c 4 x 4 + c e x 6 + ... + c 2m x 2m + 0{x 2m+1 ) 

as F xx (0) = 0. We assume that F has a non-zero term in its Taylor expansion and omit the case when 
F is C°° and has vanishing Taylor expansion at 0, i.e. is "completely flat". To find the scaling law 
for the time a trajectory spends in [— 1, 1] we have to evaluate the integral: 

dx ~ / ^-dx for r — > 



-l 



r + c A x A + c 6 x 6 + ... + c 2m x 2m + 0{x 2m+l ) J_x r + c 2 kX 2k 



where 2k is the smallest index with nonzero Taylor coefficient. Since the asymptotic behaviour of the 
time spend in [— 1, 1] is independent of the interval of non-zero length centered at we can find the 
scaling law by solving the integration problem: 



SZvTx^ dx forfce { 2 ' 3 ' 4 "--} 



To evaluate the last integral we either use substituion twice or use contour integration over the contour 
in C given by the semi-circle with radius a and the interval [—a, a]. In any case, we get 

1 2fe— 1 | 

Trrdx ~ r a* as r — ► + 



The discussion can be summarized in the following result: 

Proposition 2.2. Consider the ODE 

x = r + F(x) with F(0) = = F x (0, 0) = F xx (0) and F(x) > for x^O (8) 

isswrae i/iai F is real analytic. Define m to be the smallest exponent with nonzero coefficient in the 
Taylor expansion for F at 0, then fi§|) has a saddle-node bifurcation at (r, x) = (0, 0) with scaling law 
given by 

m — 1 

2.3 Further Remarks 

Although degenerate saddle-nodes are not "generic" in the space of C°° 1-parameter vector fields 
x = f(r,x) with /(0,0) = = /a (0,0) they still might be observed in practical applications, e.g. due 
to symmetry inherent in the system or given nonlinear parameter dependencies. 



Propositions 12. II and 12. 21 show that we should expect various power laws as scalings near a saddle- 
node bifurcation for a smooth vector field. The natural question arises what happens if we drop the 
smoothness requirements of our vector field. In particular we consider the case when the vector field 
is only continuous at (r, x) = (0, 0). 

3 Non-Smooth Saddle-Nodes 

We define a saddle-node bifurcation for C° vector fields as a topological change of the phase portrait 
homeomorphic to the case of a smooth saddle-node bifurcation. We remark that the analytical de- 
scription presented in Theorem 11.11 no longer applies and refer to [5] and [5] (and references therein) 
for analytic methods in the theory of non-smooth bifurcations. Suppose we drop the smoothness 
requirement on the parameter dependence and allow R(r) G C° (M) , then we can show: 

Proposition 3.1. Let a : (0,oo) — ► R + , a(r) € C°°((0, oo)) and a(r) — ► oo as r — ► + . Then there 
exists a function R{r), R G C°(M), such that the equation x = R(r) +x 2 has a saddle-node bifurcation 
at r = with no equilibria for r > and two equilibria for r < 0. The scaling law for this saddle-node 
is given by t ~ a(r) as r — > + . 

Proof. The proof is completely analogous to the proof of Proposition 1 2 . 1 1 wit h the additional observa- 
tion that R(r) = l/a(r) 2 e C°(R) in the current setting. □ 

Proposition 13.11 implies that any scaling law can occur if the parameter dependence on r is non- 
linear and only C . This is in contrast to the fact that in the smooth case we expect power laws. For 
the rest of this paper we shall focus on the equation: 

x = r + F(x) for x G K and r G K 



with F E C°°(R — {0}) and F E C°(R). We shall show that this case precisely gives an interesting 
"intermediate" behaviour between the cases considered so far. We introduce the notation: 



r + f- (x) for x < 
r + /+ (x) for x > 



(9) 



with / ± (0) = 0, f±(x) > for x ^ 0, /_ E C°°((-oo ! 0)) and /+ E C°°((0,oo)). Observe that © has 
a saddle-node bifurcation at (r, x) — (0,0). If we use the same techniques to investigate the scaling 
law for r — * + as in ©-(O we obtain two integrals: 



A_(r) 



•/-(*) 



-dx 



and 



A+(r) 



•/+(*) 



-dx 



Observe that if /-(a;) <C /+(x) for x — > then A_(r) <C ^4+(r) for r — > + . So we can assume 
that f—(x) ~ f+( x ) as a; — > 0. In particular we assume without loss of generality (with respect to 
investigating the scaling law) that f+(x) is given and f-(x) is its even extension to (— oo, 0). For 
simplicity of notation we shall simply denote f+(x) = /(x) for x > and /-(x) = /(x) for x < 0, 
drop the subscripts ±. 

3.1 Three Examples 

We compare three key examples of ODEs x = r + fi(x) (i = 1, 2, 3) leading to further investigation. 
They are given by 



fi(x) 



h{x) = x / 3 (x) = 



for x > 



(10) 



where we use even extensions to define the functions ft on K. The examples are illustrated in Figure 

m 

r+x for 3 different r 
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Figure 1: /(x, r) = r + x a for r = —5,0, \ and a = |, 1,2 
Integrals giving the scaling laws are: 



^(r) 



1 



-<ix 



/o r + /,(x) 
The integrals can be evaluated explicitly and the results are: 

Ai(r) = 2 + 2rln(r) - 2rln(l + r) A 2 (r) =/n ( 1 + 



forz= 1,2,3 



tan 



^a(r) 



-(* 



Table 1: Scaling laws for three examples 



Function 



Scaling Law 



Order of the Scaling 



/i(aO 
h{x) 
h{x) 



2 

ln(r) 

1 



constant 

logarithmic 

square-root 



We consider the limit r — > + to get the scaling laws. The results are summarized in Table [TJ 

Therefore we should pose the question, what exponents a for the family f a (r,x) = r + x a have 
scaling laws 0(1). 

3.2 "Quantitative Bifurcation" 

Theorem 3.2. For each r, define f a (r,x) as the even extension of 

x i— » r + x a x > 

Consider the two-parameter family of differential equations given by 

x = f a (r, x) for i€l and reR 

then we find the scaling law for a £ (0, 1) to be given by t = 0(1), whereas t — > oo for a £ [1, oo), as 
r^0+. 

Proof. Let r n > be a sequence such that r„ > r n+ \ for all n and r n — > as n — > oo. Then consider 



lim 



-da; 



(11) 



'n 



Notice that g n (x) 



is a montonically increasing sequence of functions on [0, 1]. Therefore we 



can apply the monotone convergence theorem in equation (JTTJ) to obtain 

'1 1 /■! 



lim 



dx = / lim 

+ X a J n-»oo r„ + X c 



1 1 

dx = I — dx 
x a 



10 In -T * Jo '"^inTl Jo 

where the last integral converges for a £ (0, 1) and diverges for a > 1. □ 

The main observation of Theorem [373] is that there is a major quantitative change in the behaviour 
of the solutions at a = 1 . In analogy with the classical terminology we refer to a — 1 as a "quantitative 
bifurcation" . Notice that a — 1 lies on the boundary of C° and C 1 functions spaces for the family 
f a (r,x). 



It is an easy extension of Proposition 12.21 that for a > 1 the scaling laws follow powers of r, i.e. 
t ~ r- 7 with j £ (0, 1). We have seen in Table Q] that for a — 1 the scaling is logarithmic. Hence the 
situation observed near a = 1, providing a transition from a power law to a constant via a logarithmic 
"break-point" , resembles the situation of classical qualitative bifurcation theory in the context of a 
quantitative feature of a dynamical system. 





Figure 2: Illustration of the simple pendulum with a barrier at ^ at distance 1 from the center. 

4 A Model Problem 

Consider a pendulum as shown in Figure [2j 

We use the following notation: I is the length of the pendulum, m is its mass, g is acceleration due 
to gravity, v is a viscous damping coefhcent, ft denotes constant forcing, 9 is the angle and / denotes 
the moment of inertia of the pendulum. Then Newton's law gives the equation of motion as: 

16" + mgl sind = v(u> - 6') 

If we consider the overdamped limit of very large damping we can neglect the term 19" and consider 
the classical nonuniform osciallator 

0i =UJ - r ^-sm6 (12) 

v 

For simplicity let us set L = mgl/v : then (TT2"|) reads: 

9' =lj- Lsin8 

Usually one assumes that the length of the pendulum is fixed and the bifurcation parameter is the 
applied torque. We generalize this approach and assume that we can vary the length L depending 
on the angle of the pendulum, i.e. L = L{9) with L : S 1 — ► R + . Note that we regard S 1 as [— ir,ir] 
with endpoints identified as indicated in Figure [2] If we assume that we want to vary the length 
symetrically on S 1 there will be breakpoints for L{9) at 9 = w, f ,0,— ?. Furthermore let us assume 
that the length of the pendulum is limited at one of the breakpoints, say L(tt/2) = 1 (see Figure [5]). 
Let us consider the family of functions: 

C -l + (-2/n(9 + n/2)) a for0e[-7r,-§) 

I -l + (2M9 + n/2)r for ^e [-f,0) 

qW | 1 - (-2/tt(0 - 7r/2)) a for0e[O,§) 

[ l-(2/7r(6»-7r/2)) a for €[§,*] 

for a > 0. Illustrations for a = i, 1, 2 are given in Figure [3] 



square-root distorted wave 



triangular wave 



parabolically distorted wave 




Figure 3: f(0, 1) = 1 + F a {6) for a = i, 1,2 



Now we consider the equation modelling the elongation L as: 



L(9) 



Fa(9) 



sine 



Note that this is a priori not defined for a < 1 at giving infinite length L. We can simply truncate 
L(9) if we want to construct the experiment but we remark that the form of L{9) near is not relevant 
for the following discussion. It is crucial to notice that L(tt/2) = 1 independent of a, i.e. we pass the 
bottleneck with the same length. We are left with the equation 



F a (9) for 9 e [-7T, 7r] and a, uj e 



(13) 



As in the usual sinusoidal case 9 = w — sin 6* we have a saddle-node bifurcation at (u>, 9) = (1, ? ) for all 
a > in equation (JT^J) (see also Figure [3]) ■ From Theorem 13.21 we see that there exists a quantitative 
bifurcation for a = 1. In particular, a € (0,1) gives scaling laws t ~ 0(1) as w — > 1 + . For the 
triangular wave a — I, we get a logarithmic scaling and for a > 1 we obtain power laws £ ~ (w — 1) J 

for i e (o,i). 



This means that upon tuning the parameter a in the family F a and setting uj small and positive 
we can switch the behaviour of the bottleneck. Hence equation fT3")) describes how different strategies 
of varying L affect the motion. The case < a < 1 corresponds to reducing a long pendulum to unit 
length for 9 — ► (§)~, the case a = 1 corresponds to a very small variation upon approaching J and 
a > 1 means that we start with a very short pendulum with a rapid increase near 7r/2 to reach unit 
length. 



5 Conclusions 

We have investigated scaling laws for saddle-node bifurcation in 1-dimensional dynamics. The generic 
smooth 1-parameter vector field exhibits a square-root scaling law. Dropping the assumptions of 
genericity we have seen that the nonlinearities of the parameter or the phase-space variable exhibit 
various power laws. Furthermore we have investigated the case of C° vector fields and demonstrated 
that parameter linearities lead to arbitray scaling laws. For phase space nonlinearities we found a 
quantitative bifurcation for a 2-parameter family of C° vector fields. 



In particular, the theory presented for C° vector fields can clearly be extended to more than 1 
dimension, but in contrast to the smooth cases we do not have tools like center manifolds immediately 



available. As an example it is easy to construct an n-dimcnsional system, which has locally at least 
two directions near a saddle-node with different C° vector fields in each direction. This substanially 
complicates the analysis for the natural extension to more dimensions. 

We have also demonstrated that a very simple pendulum equation can exhibit a quantitative bi- 
furcation. Many other applications of saddle-node bifurcations occuring in n-dimensional C° systems 
could clearly exhibit the phenomenon. 

Furthermore we mention that the different scaling laws found for C° vector fields can be directly 
related to the same scaling laws found in C° maps. C° maps are of particular relevance in the 
analysis of discontinuity-induced bifurcations and their associated return maps (Poicare discontinuity 
map [PDM] and zero time discontinuity map [ZDM]) [2]. These return maps have different types of 
singularities, among them we can find piecewise-linear maps and maps with square-root singularities. 
Hence the methods we presented in this paper are likely to be very useful in the analysis of scaling 
laws for discontinuity-induced bifurcations. 
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